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We study the low-temperature behavior of a magnetic impurity which is weakly coupled to 
correlated conduction electrons. To account for conduction electron interactions a diagrammatic 
approach in the frame of the 1/N expansion is developed. The method allows us to study vari- 
ous consequences of the conduction electron correlations for the ground state and the low-energy 
excitations. We analyse the characteristic energy scale in the limit of weak conduction electron inter- 
actions. Results are reported for static properties (impurity valence, charge susceptibility, magnetic 
susceptibility, and specific heat) in the low-temperature limit. 



I. INTRODUCTION 

Metals with strongly correlated electrons exhibit 
highly complex phase diagrams at low temperatures re- 
flecting a rich variety of possible ground states. Promi- 
nent examples are the well-known normal Fermi liquid 
state as well as magnetically ordered, superconducting 
and insulating phases which may coexist or compete 
within the same material. A key to a quantitative under- 
standing of the unusual phases is therefore a quantitative 
description of electronic correlations and their observable 
consequences. 

The present paper focusses on (dilute) magnetic alloys 
with correlated conduction electrons, i.e., we consider 
host metals with correlated conduction electrons contain- 
ing a small amount of magnetic ions. We investigate the 
question how conduction electron correlations affect the 
formation of a nonmagnetic Fermi liquid ground state 
commonly referred to as Kondo effect. The latter has 
been known to be the source of many anomalous prop- 
erties in magnetic alloys with nonintcracting conduction 
electrons. In addition to its relevance in magnetic alloys 
the Kondo effect is becoming important in the study of 
interacting mesoscopic systems. Theoretical techniques 
which provide a detailed quantitave understanding of the 
physical properties of these systems are hence highly de- 
sirable. To leading order in the low impurity concentra- 
tion the electronic properties of dilute magnetic alloys 
can be calculated in two steps. First one has to deter- 
mine the elctronic properties of the host which will not be 
significantly affected by the addition of a small amount 
of impurities. In the second step the contribution of the 
magnetic ions has to be calculated. 

For a metal with uncorrelated conduction electrons the 
first part of the problem is solved by standard meth- 
ods of electronic structure calculation. The theory for 
the second step is well established 1-5 . The theoretical 
techniques available include exact solutions for equilib- 
rium properties as well as approximate methods for dy- 



namic properties. Of particular importance in this con- 
text is the diagrammatic approach based upon the large- 
degeneracy expansion. This scheme can be generalized 
to the treat non-equilibrium properties which makes it a 
very flexible method. 

The central goal of the present paper to extend the 
large-degeneracy expansion for the normal-state proper- 
ties of dilute magnetic alloys to the case of host metals 
with correlated conduction electrons. In this case, the 
first step, i. e., the treatment of the host is a highly 
non-trivial problem which has not yet been solved. Par- 
tial answers, however, can be found in limiting cases. 
Although the diagrammatic approach developped in the 
present paper is valid for a general conduction elec- 
tron interaction ( CEI ), we focus on systems where 
the ground state and the low-energy excitations of the 
interacting conduction electrons smoothly evolve from 
those of the non-interacting reference system. This is 
in marked contrast to the specific behavior encountered 
in one-dimensional systems. Theoretical studies which 
have been performed for various models including both 
impurity spins 6-8 and Anderson impurities 9 11 coupled 
to Luttinger liquids predict rich phase diagrams. Adopt- 
ing well-established models for the electronic properties 
of the host, we calculate the evolution of the charac- 
teristic energy Tr of the low-lying magnetic excitations 
with the conduction electron repulsion. In the case of 
uncorrelated conduction electrons the latter is usually 
much smaller than the typical energy scale of the con- 
duction electrons set by the band width D and depends 
exponentially on the inverse coupling between the local- 
ized electron and the extended conduction states. This 
fact is a direct consequence of the Fermi liquid ground 
state realized in normal metals. The diagrammatic ap- 
proach allows us to explicitly and quantitatively study 
how the different consequences of electronic correlations 
(mass renormalization, effective interactions etc) affect 
the Kondo effect. 

The main scope of this paper is to analyze how CEI 
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influence the contribution of magnetic impurities to mea- 
surable properties in general and its scaling properties 
in particular. We calculate thermodynamic properties 
(impurity valence, charge susceptibiliy, magnetic suscep- 
tibility and specific heat) in the low-teperature limit to 
leading order in the inverse degeneracy. 

Recent calculations for a magnetic impurity in a metal 
with interacting conduction electrons 12 ' 13 adopted either 
the DMFT approach 14 or the NRG but for a very special 
model 15 . The model calculations mentioned above pre- 
dict non-trivial variation with the Coulomb repulsion of 
the characteristic temperature Tk- 

Generally, the modifications introduced by the conduc- 
tion electron interactions (CEI ) into the low-energy ex- 
citations arise from the subtle interplay of three different 
types of influences. First, the density of conduction states 
at the Fermi level is changed. Second, the probability 
for virtual transitions between impurity and conduction 
states are reduced by the on-site Coulomb interaction U. 
Third, the effective spin coupling between the conduc- 
tion and impurity electrons is enhanced by the increased 
number of uncompensated spins in the correlated con- 
duction electron system. Considering these facts, it is 
not surprising that model studies accounting only for se- 
lected aspects arrive at rather controversial conclusions 
concerning the Kondo effect in metals with correlated 
electrons 16-18 . The Kondo spin model generalised to 
the case of the interacting conduction-electron host was 
discussed in 16 and it was shown there that two-particle 
Green's functions of host electrons ( vertex corrections 
) are an essential component of the theory which leads 
to an enhancement of the exponential Kondo scale for a 
weak CEI. This enhancement may be traced to the third 
type of effects caused by the CEI. The ground state en- 
ergy of the Anderson impurity for weak CEI was consid- 
ered in the frame of 1/N expansion 17 . The same enhance- 
ment of the exponential Kondo scale, formally due to 
the renormalization increase of the hybridization width 
A, appears in this work 17 . In contrast to the above- 
mentioned findings a decrease of Tk due to the CEI in 
the Hubbard model was reported in the paper 18 . This 
decrease is a consequence of the change in the single elec- 
tron properties of conduction electrons caused by the in- 
teraction U ( including the change of the chemical poten- 
tial as the function of U ) . The vertex corrections influ- 
ence which renormalizes both parameters of the Ander- 
son impurity model 19 , e/ and A are not considered in 18 . 
At this point, we should like to mention that the role 
of the Coulomb interaction between the magnetic impu- 
rity electron and conduction electrons, Uj c . was broadly 
discussed. We do not discuss here the Coulomb interac- 
tion between localized and conduction electrons which is 
considered in its various aspects in 20-24 . It was shown 
that its effect at Uf c <C Uf may be fully absorbed by 
the renormalization of two parameters of the Anderson 
impurity Hamiltonian: the impurity electron energy level 
Ef and the hybridization width A. In the following we 
assume that the on-site impurity electron Coulomb re- 



pulsion Uf is very large, Uf — > oo, and we do not take 
into account explicitly the Uf c interaction. 

The paper is organized as follows: In Section II we 
begin with a discussion of the Hamiltonian for an An- 
derson impurity embedded in a metallic host with cor- 
related conduction electrons and the extension of the 
standard selfconsistcnt large-degeneracy approximation 
to the case of CEI. Both the interaction-induced changes 
in the single-electron spectral function of interacting con- 
duction electrons and their vertex function are included. 
In Section III expressions for configurational selfenergies 
together with the NCA integral equations are formulated 
for a general case of CEI. The expressions are evaluated 
in Section IV for a model where the Coulomb vertex 
function is only weakly frequency-dependent. Thermod- 
inamic properties at zero temperature are presented in 
Section V and Section VI contains discussions and sum- 
mary. Technical details related to the explicit evaluation 
of diagrams and discussions of the fourth-order hybridis- 
ation coupling are in the appendices. Some of the results 
appeared in the short unpublished preprint 19 . 

II. MODEL AND CALCULATION AL SCHEME 

We adopt a generalized Anderson model for a mag- 
netic impurity coupled to interacting conduction elec- 
trons. The resulting Hamiltonian reads 

H = Hband ~t~ Himp ~\~ Hmix (1) 

where the three components describe the conduction elec- 
trons, the f states and a hybridization or mixing interac- 
tion between the two 

Hband = e k C l a C ka + H CEI 

ka 

Hi m p — ^ ' £/^/m H ^ ' Tlf m nf m i 

H mix = Y, (VmAk)fLc^ + h.c.) . (2) 

k,m,(T 

The creation (annihilation) operators for conduction elec- 
trons with momentum k, band energy eg and spin a are 
denoted by (c£ CT ). Throughout this paper, all energies 
are measured relative to the Fermi level. The conduction 
states are assumed to be orbitally non-degenerate. Their 
interaction is accounted for by 

Hcei = Yl u $+ <f> *' - <?; *) 

k,k',q a,a' 

xct ct c,-, ,Cr ■ (3) 

k+qa k'-qa 1 k <7 ka \ > 

where L is the number of lattice sites. In the present 
paper, we approximate Hcei by a Hubbard-type inter- 
action 
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U(k + q,k' -q;k',k) -» U 



(4) 



where U denotes the local Coulomb repulsion between 
two conductionclectrons at the same lattice site. Another 
important example which shall be studied in a forth- 
coming paper are Fermi liquid systems where the CEI 
renormalizes the quasiparticle dispersion eg and also in- 
troduces a 'residual' interaction among them. 

The fJnifm) are the creation (annihilation) operators 
for /-electrons at the impurity site. They are charac- 
terized by the total angular momentum J and a quan- 
tum number m which denotes the different states m = 
1, . . . , N within the JV-fold degenerate ground state mul- 
tiplet with orbital energy e/. The Coulomb repulsion 
Uf between two /-electrons at the impurity site is as- 
sumed to be much larger than the other energy scales 
and therefore we may let Uf — ► oo. For simplicity we do 
not include here excited multiplet states, ignore crystal 
electric field splittings and assume that the impurity has 
only one electron ( hole ) in its magnetic configuration. 
We account for the large Coulomb interaction among the 
f-clcctrons Uf — > oo by restricting the Hilbcrt space, i. e., 
by removing all states in which the f occupancy exceeds 
unity. 

The mixing between the two subsystems is conve- 
niently characterized by the "hybridization width" 25 
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We are mainly interested in the regime |e/| 3> A TO (0) 
which is usually referred as "local moment regime". 
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FIG. 1. Self-consistent f configuration self-energies and 
contributions to the vertex. The solid, dashed and wavy 
lines represent the dressed propagators for conduction elec- 
trons, occupied and empty f states. The open circle denotes 
the bare hybridization V while open and filled squares are 
the bare on-site Coulomb repulsion and the two-particle ver- 
tex r* C ^(l, 2; 3, 4), respectively, (a) Empty state self-energy 
£(i(i^ m ). (fi) Occupied state self-energy £o-(icj n ). (c) Contri- 
bution to the effective hybridization vertex, (d) Lowest order 
correction. 



The central goal is to calculate the impurity contribu- 
tion to the low-energy properties of the dilute magnetic 
alloy. The latter are given in terms of the Green's func- 
tions for the empty f-state |0 > (4/°- configuration) and 
the occupied / states \m > (4/ 1 - configuration) denoted 
by Gq(z) and G m (z), respectively, 



G (z) 



1 



z — 



£o(z) 



G m (z) 



z-ef- T, m (z) 
(6) 



They are coupled through the configurational selfener- 
gies Tiq(z) and and X m (z) for which we derive expressions 
proceeding in close analogy to the case of non-interacting 
conduction electrons. 

The electronic properties of the metallic host are not 
affected by the presence of a small amount of magnetic 
impurities. To leading order in the small concentration 
they are characterized by the 1- and 2-particle Green's 
function describing the single-particle excitations and the 
two-particle correlations of the interacting conduction 
electrons, respectively. 

We assume that the single-electron Green's function 



G a (k, z) — 



1 



Z) 



(7) 



as well as the conduction electron selfenergy z) do 

not explicitly depend upon the wave vector k but vary 
with k mainly through the bare band energy, i.e. 



G a (k, z) — G CT (eg, z) 



£o-(fc, z) = £ CT (e£, z) 



(8) 



This condition is always satisfied in the DMFT 
approach 14 where the dominant many-body effects are 
included in a local selfenergy. As a consequence, also 
the general n-particle Green's functions of the conduction 
electron system depend upon the wave vectors through 
the corresponding band energies. 

The configurational self-energies So and S m are de- 
rived by means of a perturbation expansion in terms 
of Green's functions for the f-configurations, 1- and 2- 
particle Green's functions for the interacting conduction 
electrons as well as (bare) hybridization vertices. The 
rules for constructing and evaluating the empty- and 
occupied-state selfenergies in the restricted Hilbert space 
of the infinite- U- Anderson model are summarized in 26 . 
Typical contributions to the f-configurational selfenergies 
are displayed in Figure 1 . These include the non-crossing 
diagrams Figure 1(a) und (b) where the conduction elec- 
tron interactions enter through the fully renormalized 
conduction electron propagator. The diagrams 1(c) and 
(d) describe vertex corrections. We shall show below 
that under the assumptions Eq. (8) the infinite-order 
summation of these diagrams based on the self-consistent 
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approximation for the empty- and occupied-f-state prop- 
agators can be considered as the leading order contribu- 
tion in the inverse degeneracy 1/N. The remaining part 
of the present section is devoted to the justification of 
this conjecture. 

We start by briefly summarizing the basic facts on 
which the large-degeneracy expansion is based in the case 
of non-interacting conduction electrons. The classifica- 
tion scheme exploits the fact that the bare conduction 
electron propagator l/(z — eg) depends upon the wave 

vector k through the bare band energy. As a conse- 
quence, the summations over internal fc-vectors can be 
decomposed into integrals over the (bare) band energies 
e and averages over constant energy surfaces according 
to 

lE"-^ j ^\H^-^)--- = JdeN(e)(...) e 

k k 

(9) 

Here N(e) is the density of bare band energies. The k- 
averages (. . .) e which contain only combinations of the 
hybridization matrix elements 



< V^ a (k)V m , a (k) > £ = 



(10) 



1 



1 



V(f) | Z ^C CT (fc)Wfc)^-eg) 

k<7 



provide the m-selection rule 27 which simplifies the struc- 
ture of the selfenergy contributions and ultimately allows 
for a classification with respect to the small parameter 
1/N. 

^From the preceding discussion it is apparent that 
the assumption Eq. (8) guarantees the validity of the 
1 /N classification scheme for all diagrams where the con- 
duction electron properties enter via the single-particle 
Green's function. Within this subclass the contributions 
displayed in Figure 1(a) and (b) ( without vertex cor- 
rections ) are the leading ones with respect to the small 
parameter 1/N. 

To assess the validity of the 1/N expansion is more sub- 
tle for the diagrams containing the two-particle Green's 
function. Here the simplifying assumptions Eq. (8) im- 
ply ( see Section III ) that the hybridization matrix ele- 
ments enter diagrams in Figure 1( a,c ) and ( b,c ) in the 
combination 



W 4 )(£i<7i, E 2 <72; E 3 a 3 , £404) = — 

ki ,k2 : kz,k& 

S(E 1 - e(h))6(E 2 - e(k 2 ))6(E 3 - e(k 3 ))5(E 4 - e(fc 4 )) 

V ™°1 (kl)V m 'a 2 (fe)KL 3 (>3)*C'<74 (M 



S* (k 3 + k 4 - h - k 2 ) . 

where the Lauc function 
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S*(ki + k 2 -k 3 -k 4 ) = i ^2 CX P {* (ki + fa-fo-kij j 

(12) 



Rn 



accounts for momentum conservation up to a reciprocal 
lattice vector. 

In the Appendix A we present a detailed model cal- 
culation for a rare earth impurity hybridized with tight- 
binding s-band states. The results show that the new 
contributions to £0 and T, m are 0(1) and 0(1/N), i. e. , 
of the same order of magnitude with respect to 1/N as 
their NCA counterparts. It is interesting to note that 
the dominant contributions to are non-local com- 
ing from the coupling of the /-states to the conduction 
electrons at the neighboring sites. 

To summarize, the configurational selfenergies dis- 
played in Figure 1 provide a consistent extension of 
the well-known selfconsistent large-degeneracy expansion 
to the case of interacting conduction electrons. Ac- 
tual calculations, however, require the fully rcnormalizcd 
conduction-electron propagator as well as the Coulomb 
vertex. Since this problem still remains unsolved for the 
Hubbard model 28,29 we have to adopt approximate ex- 
pressions derived either from phenomenological consid- 
erations or from partial resummation of selected classes 
of diagrams. 

General qualitative results can be derived in limiting 
cases. Prominent among them is the case where the 
Coulomb vertex can be considered as a static quantity 
which includes the limit of weakly interacting conduction 
electrons as well as the Fermi liquid case. 



III. CONFIGURATIONAL SELFENERGIES 

For non-interacting conduction electrons, the self- 
consistent solution 30 ' 26 has three characteristic features: 
The occupied f-spectrum shifts to peak at a value e/ — 
e/+3?S m (e/), the dominant contribution to the level shift 
coming from the continuum of charge fluctuations. The 
resonance in the occupied f-spectral function acquires a 
small width. Finally, the empty state spectral function 
exhibits a pronounced structure at ljq — If — Tk which 
develops with decreasing temperature and which sets the 
scale for the low-temperature behavior. This feature is 
the direct manifestation of the Kondo effect reflecting 
the admixture of /°-contributions to the ground state 
and the low-energy excitations. 

In this paper, we study the influence of the CEI on this 
non-perturbative feature. Of particular interest are the 
position of the resonance energy luq relative to the energy 
if of the 4/ 1 configuration as well as the weight of the 
resonance. 

Let us first neglect vertex corrections and focus on 
the modifications introduced by the CEI into the single- 
particle excitations of the conduction electron system. 
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They are accounted for by inserting the full conduction 
electron propagator for interacting electrons G a {k,iuj n ) 
from Eq. (7) into the configurational selfenergies Figure 
l(a,b). 

The selfenergy of the occupied f-level 
S^M =^£W*)/ dE n f (-E)A„(k,E) 

ka 

xG (u-E)V* a (k) 
= lJdE A° m (E)n f (-E)G (uj E) (13) 

is diagonal in m as shown in the previous section. Here 
n f (E) denotes the Fermi function. The properties of the 
metallic host are reflected in the energy-dependent hy- 
bridization strength 

A ln } (E) = ^V m(T (k)A„(k,E)V^(k) (14) 
No- 
where the conduction electron spectral function 
A„{k,E) = -~±ImG a {k,E) depends upon the wave vec- 
tor fc mainly through the band energy eg ( see Eq. (8) ). 
The selfenergy of the empty state is treated in the same 
manner so the corresponding selfenergy expressions re- 
duce to 

4 0) (a,) = - J2 [dE A^(E)n f (E)G m ( W + E) 

H = I J dE A<£\E)n f (-E)G ( w - E) (15) 

in close analogy to the case of non- interacting electrons 30 . 

The selfconsistency equations Eq. (15) were solved 
numerically 31 for various well-established approxima- 
tions to the spectral function of interacting conduction 
electrons such as the Hubbard III model 32 and the Roth 
approximation 33 . The general results can be summa- 
rized as follows: It is obvious that for (weakly) inter- 
acting conduction electrons the dominant effect of hy- 
bridization on the 4/ 1 configurational spectrum is a shift 
ef(U) — e/ ~ Jt£ m (e/) of the resonance energy which is 
renormalizcd by the Coulomb repulsion U and its influ- 
ence on the charge fluctuations. The quantity of interest, 
however, is the empty-state selfenergy and its variation 
with energy in the vicinity of e/ which can be deduced 



from rather simple considerations assuming that the CEI 
do not introduce anomalies into the conduction electron 
spectral function on the energy scale defined by the char- 
acteristic temperature Tk- The smooth variation with 
energy of J2k A<r{k, E) implies that in the metallic state 

the basic analytic structure of y!^\uj) is not altered as 
compared to the case of non-interacting conduction elec- 
trons, the characteristic feature being a logarithmic vari- 
ation in the vicinity of the f-energy e/. The prefactor, 
however, is proportional to the interaction-renormalized 
density of states at the Fermi level N(ef). The low- 
energy scale, Tk, i.e. the distance between the pole in the 
empty f-state Green's function and the 4/ 1 peak depends 
on the renormalized parameters in the usual exponential 
way. Especially the above is clear for the case when the 
CEI leads to the spectral function of the quasipartical 
type, A a (k, E) = S(E — eg) with a new dispersion eg. 

Would the single-electron contribution presents the 
whole story the CEI-case would be relatively simple. The 
Coulomb interaction induces vertex corrections which are 
of the same order in the inverse degeneracy l/N as the 
preceding single electron contributions. They are an im- 
portant ingredient of the theory and must be included in 
the discussion 16,19 . 

The explicit calculation requires the full Coulomb ver- 
tex T^ 7 ) of the conduction electrons as input which must 
be determined consistently with the conduction-electron 
self-energy. We evaluate the vertex corrections by ana- 
lytic continuation from the Matsubara frequencies insert- 
ing the spectral representation 

G a (k,iw m )= [ dE MhE l (16) 
J iuj m - E 

for the conduction electron propagators and following the 
rules specified in 26 . The projection onto the relevant 
physical subspace is performed implicitly in the summa- 
tion over the Matsubara frequencies where we retain only 
the contributions from the poles in the conduction elec- 
tron propagators. The empty state self-energy, (w), 
(see Figure (1, a, c))can be written as 

sf J H=^E fdEAW(E,u)n f (E)G m (u> + E) (17) 

m 

where the Coulomb contribution to the hybridization 
strength is given by 



A^l(E,co) =7 ri^A CT (fc,£)-l^ J <Kl n f (-E + n)A ai (k u E-n)G (uj + n) 

ka fei(Ti 

EE E / dE ' Gm> (^ + E ' + °) 

{A a2 {k 2 ,E' + Q)n f {E' + n)G a , {%' ',£?') + A a ,(k' ',E')n f (E')G a2 (fe,£7' + n)} 
(VU (h)V m >a 2 {k2)V: ia {k)V mla ,{k')*Yf^ aa , (hE - n,k 2 E' + 0; kE, k'Ef) + c.c) 
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S* (k! + k 2 -k- fc') 



(18) 



Here and elsewhere T < ^l 2 - (73C r i (1, 2; 3, 4) is the for the hybridization strength entering the occupied-f- 
Coulomb vertex corrections with indices 1,2 for in- and states selfenergies, ( see Figure (1, b, c ) 
3,4 for outgoing particles. A similar expression is found 

s £° M = I J dE A^ m (E,uj)n f (-E)G (uj - E) (19) 

with 

AM m (E,u) ^E^'^EE j dn n f (E-Q)A ai (k u E-n)G m , (w - SI) 

ka kicn m ' 

EE / dE ' G ° ( w - E ' - n ^ 

k 2 a&' a' 

[A a2 {k 2 , E' + Sl)n f {-E' - Sl)G a , (k 1 , £') + A a ,{k', E')n f (-E')G„ 2 (k 2 ,E' + fi) } 
(v ma {k)V m , a ,(k')V^ ai {k)V^ a2 (fc 2 )r^ ;CTiCT2 (kE, k'E'; k x E - SI, k 2 E' + Sl)+ c.c) 
S* ^+k 2 -k- fc') (20) 



Note that the self-consistency equations, Eq. (15) gen- 
eralized by including the vertex correction contributions 
from Eqs. (18, 20) in the integrands of Eq. (15) read 

SoH = ^E/ dE(KHE) + 

m 

A U i(E,cj))n f (E)G m {w + E) 

E m {u) = ±J dE (A^(E) + 

AM m (E,Lo))n f (-E)G Q (u-E) (21) 

Eqs. (18, 20) are general in the sense that they do not 
assume any specific form of conduction electrons spectral 
functions, vertex corrections, etc.. In the case when k- 
dependences in conduction electron propagators enter as 
in Eq. (8) only via the conduction eletrons dispersion 
e(k) Eq. (18) may be casted in the form 

A U m(E, w)) = J dojidoj 2 <kjzduj i ^ 

C\a 2 a 3 a^ 

V^;\criUJi, a 2 u> 2 , a 3 uj 3 , OiUJi) 
F m (E,u>,aiUi,a 2 w 2 ,a3U)3,<T4U>4) (22) 

Here Vm^ denotes the 4-order hybridization 
coupling given explicitly in Eq. (11) while 

F m (-E/ 1 , to , Cito o 2 co 2 , , (J4CJ4) is readily obtained by 
using Eq. (8). Eq. (22) is simplified enormously in the 
case when conduction electrons spectral functions may 
be approximated by the quasiparticlc spectra . A similar 
coupling Vm^ may be introduced for the Amm(E, oj). For 



a model calculation of the 4-ordcr hybridization coupling 
Vm^ see Appendix A. 

To summarize, we generalized the selfconsistent large 
degeneracy expansion to the case of correlated conduc- 
tion electrons. The modifications due to the interaction 
enter via the spectral function of the conduction electrons 
as well as an effective renormalized hybridization vertex. 
The explicit evaluation hence requires these quantities 
for a system of interacting conduction electrons. In the 
subsequent sections, we shall consider the influence of 
the Coulomb repulsion U on the effective hybridization 
strengths which depend upon both E and u>. In particu- 
lar, we shall discuss the analytic structure of the selfener- 
gies for weakly interacting electrons and discuss the mod- 
ifications in observable properties in the low-temperature 
limit T -> 0. 



IV. WEAK CONDUCTION-ELECTRON 
INTERACTION 

As a first example, we consider the limit of weakly 
interacting electrons, i. e. , we assume the Coulomb re- 
pulsion U to be much smaller than the bandwidth D. To 
leading order in the small ratio U/D « 1 we can ne- 
glect changes in the spectral function of the conduction 
electrons which we assume to be given by 

A a {k,E) -^S(E-e(k)) . (23) 

The central focus of the present paper is the lowest pole 
wq of the cmpty-f state Green's function 
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- 4 0) M - 4 U) M = (24) r%>. nn (kE, k'E'; k\E fi, k 2 E' + ft) 



(26) 



and its variation with the Coulomb repulsion U. 

The Coulomb repulsion contributes to the config- 
urational sclfenergies via the effective hybridization 
strengths Eq. (18) and Eq. (20) where the full Coulomb 
vertex is replaced by the bare local non-retarded repul- 
sion ( see Figure 1, d ) 



for the case of Am m accordingly. 

We elaborate on the selfenergies expressions, Eqs. (17, 
19), for the case of an orbitally non-degenerate Ander- 
son model. In this case the hybridisation matrix element 
reduces to V ma (k) = 5 ma V(k) and the occupied f-state 
propagator does not depend upon the m-indcx 34 . 

Inserting Eqs. (23, 25, 26) into the vertex corrections 
Eq. (18) and Eq. (20) correspondingly we obtain 



CT 1 CT 2 '.(7(7 



k x E - n,k 2 E' + fl; kE, k'E' 



(25) 



U 2 



( i( ry)c 



for the case of A^ and 



= ttU [ dE 1 dE 2 dE'nfi-E^Go (uj + E — E x ) 

J Ei 



2SRV ( - i HE 1 ,E 2 ,E,E') 



+ E 2 - E-E' 



and 



{G m (w + E 2 )n f {E 2 ) - G m (uj + E' + E - Ex) n f (E')} 

2^tV^(E,E',E 3 ,E 4 ) 



(27) 



A%) m (E,uj) =ttU J dE' dE 3 dE 4 n f (E 3 )G- m (uj + E 3 - E) 



irU 
{Go (w 



E 3 + E 4 -E-E' 
E 4 )n f (-E 4 ) -G (lv + E 3 -E- E')n f (-E')} 



(28) 



To derive these expressions we used Eq. (11), per- 
formed the cr-summations and the relevant integrations. 
The hybridization matrix elements V? vary smoothly 

with k , and, as a consequence, Vm\E\, E 2 , E 3 ,E 4 ) is a 
smooth function of the energies Ei i = 1, ... ,4. It can 
be approximated by (see Appendix A) 



and 



v^(e u e 2 ,e 3 ,e 4 ) 



\2irD J 



(29) 



where A is the hybridization width. In the following we 
shall adopt a flat density of states extending over the 
energy range (—D,D) and use for the vffl the last ex- 
pression in Eq. (29). 

We start by discussing the configurational sclfenergies 
for T = where the Fermi function can be replaced by 
the step function nj(x) = 0(—x) and we insert the free 
propagator for the occupied-state Green's function 35 



G m (uj) 



LU-e f 



(30) 



The selfenergies T, U ^ and can be expressed in terms 
of three integrals I oi and I mi , 1 = 1,2,3, respectively, 



-ioiM m 



uj — e f — D 



+ I 02 (uj) + I 03 (uj) 



(31) 



-J m i(w)ln 



uj 



uj- D 



I m 2{u) + I m3 {oj) 



(32) 



Further we discuss the half- filling case 37 . This particular 
choice of the band filling, however, does not affect the 
analytic behavior in the energy range of interest, i. e. , 
for uj ~ ef. 

We should like to emphasize that the integrals Ioi and 
I rn i in Eqs. (31) and (32) depend upon the full empty- 
state Green's function Gq(u>). This fact implies that the 
Coulomb contribution to the selfenergy, E^(u;), has to 
be determined selfconsistently from Eq. (31) in prin- 
ciple. In the present paper, we employ an iterative 
scheme and adopt a convenient parametrization of the 
spectral functions Aq(uj) and . Before presenting the re- 
sults, let us briefly summarize our procedure. In the 
first step, we insert the free empty-state propagator, i.e., 
A (uj) — > S(uj) into the rhs of Eq. (31). The resulting 

selfenergy T, q U \uj) yields a Green's function Gq(uj) which 
has a Kondo-type pole at ujo c (U) < ujq(U = 0) < e/ 
with rather small weight 1 — rif c (U). The index c de- 
notes the fact that only the charge fluctuation contri- 
bution was included in the selfenergy Sq^(w). In the 
next iteration, we account for the low-energy peak in 
the spectral function which we model by two (^-functions 
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4>(w) - (1 - n /c (^))5(w - woc(EO) + n/ c (^)5(w). In- 
eluding the low-energy spin fluctuations furthers shifts 
the threshold u>o(U) to lower energy, i.e. we find ujq{U) < 

UJ 0c {U) < LU (U = 0) < £/. 



Modelling the spectral function A (uj) by a combina- 
tion of <5-functions allows us to decompose the integrals 
into contributions from the charge fluctuations Ioi C (u)), 
Imic{u) ( further all paramerers and variables are in units 
of the band half-width D ) 




rl rl 2 In 

ioicM = / dx dy- — 

Jo Jo (-w+a 



-oj+ef+x+y+1 



rl In 



Io2c{u) = dx 



— LJ + X-\-l 



+x+y)(-Lu + e f +x) 

2 



Zd3 C ( w ) = / da;/ dy / <iz21n 
Jo Jo Jo 



x + y + z + 1 



x + y + z 

X [(-LJ + X + y)(-U> + €f + x) 

x (-cj + ez+x + y+z)] -1 



/mlc(w) = dx 

Jo Jo 



2 In 



-u+x+y+1 



4i2c(w) 



o Jo (-u+x)(-u}+e f +x+y) 



— UI+l + €f+X 

(-oj+x) 




) = [ dxj dy[ dz2 In 
Jo Jo Jo 



.t + y + z + 1 



.t + y + z 
x l(-LU+x){-LU + e f + x + y) 
x (-cj + .T + y+z)]~ 1 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 



FIG. 2. Contributions to the occupied-state selfenergy 
from the Coulomb correction to the effective hybridization 
vertex: Occupied-state integrals / m «c(e/) evaluated for vari- 
ous values of the f-level energy e/. Here and in all further 
Figures solid lines are guides for the eye. 

The spin fluctuations integrals are of analogical prop- 
erties but for w < uj - The infinitesimal imaginary parts 
-H0+ in denominators of the integrands in Eqs. (33-38) 
are omitted because they do not contribute for uj < e/. 
Note that for rif < 1 integrals I Mc , I nis have to be in- 
serted in Eq.(31) being multiplied by nj or 1 — n/ cor- 
respondingly. The contributions from the spin fluctua- 
tions to the Coulomb rcnormalization of the occupied- 
state vertex 7 m , s / are neglected. The Coulomb con- 
tributions to the occupied-state selfenergy vary rather 
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and from spin fluctuations integrals Ioisf and I m isf ■ The 
latter integrals are obtained from their charge fluctu- 
ations counterparts by the substitution uj — > uj — ujq. 
The charge fluctuations integrals have no singularities for 
u < £f < and it is evident just from their inspection 
that Zoic j Imic < and other integrals are positive. 



FIG. 3. Charge fluctuation contribution to empty-state 
selfenergy integrals Jo» c (w) for e/ = —2/3, A = 0.2, U = 0.1 
For uj > e/ the real part is shown. 

They give rise to a rather small shift of the effective 
f-level which can be estimated from the integrals I m ic{cf) 
displayed in Figure 2 and the real part of the selfenergy 
T,m\tf)- As we shall see below, we need not explicitly 
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account for the shift in the determination of the many- 
body low-energy scale. 

Let us now turn to the empty state selfenergy. Follow- 
ing the iterative procedure we include the charge fluctu- 
ations in the first step. The variation with energy of the 
integrals Ioi c (u>) is displayed in Figure 3. 

A detailed analysis shows that in the energy range of 
interest the integrals vary approximately like Iqi c (lo) ~ 
At In \ui - ef \ + Bi and io2cM ~ {M In \ui - e f \ + B 2 ) 2 
where A\ > A 2 . The resulting real part of the empty-f 
state self-energy varies like —(log |e/ — lo\) 2 in the vicin- 
ity of the (renormalizcd) f-lcvel. As a consequence, the 
corresponding Green's function Gq{ui) always exhibits a 
pole at ojo(U) < uio(U = 0) < ef. For not too small 
values of the hybridization width they are well described 
by the linear dependence 



to (U) ~ u (U = 0) + 



1 



(39) 



The change in the pole is seen to be proportional to the 
weight of the ^-configuration in the ground state times 
the Coulomb repulsion among the conduction electrons. 

For a first qualitative understanding of the vari- 
ation with the Coulomb repulsion of the threshold 
energy w one may use "on-shell" approximation 19 . 
Within this approximation, the empty-state selfenergy 

T,o(u>) — ^ln[e/ — cj] has the same lu- dependence 
as in the non-interacting case but with renormalizcd 
parameters: A = A (l - j^AI 01c (u> )) ;e f = e f - 

¥£r (I 02c (wo) + Io3c (w )) with luq = lu (U = 0) here. 
We see that A > A and |e/| > |e/|. If the impurity 
valence is close to integer which lead to a Kondo regime 
the renormalization of the hybridization coupling prevails 
the renormalization of the f-lcvel energy resulting in an 
effective enhancement of the Kondo energy scale. The 
simplified approach, however, cannot be used for quan- 
titative estimates. Unfortunately the variation with U 
of the corresponding Kondo-type pole is systematically 
underestimated ( Tk is overestimated ) as can be easily 
seen from the slopes 



dU 



%fr = "2^A 2 (/ 2c(Wo) + /03c(W0)) < 



The results for the Kondo pole change only slightly 
upon iteration. Inclusion of the spin fluctuation contri- 
bution to the Coulomb correction yields a rather small 
shift in the selfenergy which further stabilizes the Kondo 
ground state. This can be seen from the calculated 
variation with energy of the integrals Ioi S f(u>). The full 
selfenergy Sq( w ) is shown in Figure 4. 




> NCA+Coulomb contribution 
• Coulomb contribution 
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FIG. 4. Variation with energy of the real part of the full 
empty-state selfenergy KE^^tj) from NCA plus Coulomb 
correction to the hybridization vertex for e/ = —2/3, A = 0.2 
and U = 0.1 

The characteristic energy scale for low-energy excita- 
tions, i. e., the Kondo temperature, is now calculated 
as the difference between the ground state energy - the 
threshold u>q - and the energy of the f- level e / 



k R T, 



B 1 K — e f — ^0 



(40) 



At this point we should like to add a comment concerning 
the choice of ef. This quantity enters Eq. (40) explicitly 
as well as implicitly through u)q. If we were to account 
for the Coulomb renormalization we would have to do it 
consistently. This means we would have to consider the 
difference 

e / + Se f ~ w o(e/ + 5e f ) ~ e f - cj a + Se f (1 - n f ). The 
correction from the Coulomb contribution is hence pro- 
portional to Sef (1 — nf) which is rather small in our case 
because the shift 9?S TO (o;) is very small. 



1.5 



Charge fluctuations 
Charge + spin fluctuations 




FIG. 5. Enhancement of the Kondo temperature for 
e/ = —2/3 and A = 0.2. Charge fluctuations are dominant. 
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To conclude, we can calculate IzbTk from Eq. (40) 
using the the bare f-level energy at e/. The results dis- 
played in Figure 5 show the enhancement of the Kondo 
temperature due to the Coulomb interaction among the 
conduction electrons. 

Following up the iteration procedure the spectral func- 
tion 



A (u 



7T V 



y(0) 
^0 



Eg?)' 



(41) 



is inserted in to Eq. (31). It comes from calcultions that 
including the spin fluctuation peak at u = uj n (U) does 

not significantly alter the empty-state selfenergy £q^(o>) 
in the Kondo regime. The spin fluctuations lower the 
energy of the pole in the Green's function Go(u>) and 
therefore further stabilize the Kondo ground state. This 
can be seen from Figure 5. The data suggest the charge 
fluctuation contribution already gives a rather good esti- 
mate of the Coulomb correction to the low-temperature 
properties to leading order in the inverse degeneracy. 

The main feature of the above caculations is the factor- 
ization of the 'NCA-bubble' self-energy Sq(w) in the r.h.s. 
of Eq. (31) for the empty state self-energy T,^ (w). This 
factorisation is due to the possibility, as it was shown for 
the orbital degeneracy case in Appendix A, to neglect the 
momentum conservation in the integrals of the cou- 
pling, Eq. (11). Therefore results of this section are also 
valid for the degenerate case if one replace in Eq. (31) 
A — > (A— 1) A. So the renormalisation of the paramerers 
of the Anderson impurity Hamiltonian is selfconsistent, 
in the spirit of the NCA. 



V. THERMODYNAMIC PROPERTIES AT ZERO 
TEMPERATURE 

The results of the preceeding section allow us to assess 
the influence of the conduction electron Coulomb repul- 
sion on the thermodynamic properties of dilute magnetic 
alloys. Of particular interest are the low-temperature 
f- valence n/, the /-charge susceptibility Xo the /-spin 
susceptibility Xs an d the magnetic contribution to the 
linear coefficient of the specific heat 7. Previous calcu- 
lations based on the symmetric Anderson model yield 
a rather strong depression with U of the / magnetic 
susceptibility 13 . Data for the U-dependence of the f va- 
lence n/ and the f charge susceptibility Xc> however could 
not be obtained from these model studies since particle- 
hole symmetry pins nj to unity. For a first quantitative 
estimate we approximate the empty-state selfenergy E 

by 



J 0c 



(42) 



keeping only the charge fluctuation contribution. This 
procedure should be justified in the Kondo limit where 
the deviation from integer f valence is small. The pole 



of the corresponding Green's function, u>q, can be inter- 
preted as the ground state of the system. It yields the 
dominant low-temperature contribution to the partition 
function, and the thermodynamic properties follow by 
straightforward differentiation 5 



n f = 



du 
de f 
d 2 oj 
de 2 



lim 



7 = 



lim — m 
t->o \ T 



d 2 L0 \ 

dH 2 ) 
1 duo(T) ^ 



dT 



(43) 
(44) 
(45) 
(46) 



Finally, we shall discuss the Sommerfeld- Wilson ratio R, 
i. e. , the ratio of the zero-temperature spin susceptibility 
and the specific heat coefficient 



R 



(47) 



-7 



where fj, 2 = j (j 



l)(.9Ms) 2 . 

Our main interest is in the linear in U corrections to 
the experimental quantities. These contributions can be 
easily obtained from the linear in U corrections to the 
ground state energy as given by Eq. (39). We specify 
the interaction related enhancement/reduction in terms 
of the coefficients 



n f (U) ~ nf (l + Unf} 

Xc(U) * x< 0) (1 + u x P) 

lf{U)^lf ] (l + C/7^) 
R(U) ~ i? (0) (l + t/i? (1) ) 



(48) 



which depend upon the f- level position e/ and the 
hybridization width A and, concomitantly, on the 
Kondo energy ksT^ of the reference system with non- 
interacting conduction electrons. 

The explicit evaluation requires the generalization of 
to low but finite temperatures and to small external 
magnetic fields. The former is easily achieved by starting 
from Eq. (27) and proceeding in close analogy to the zero 
temperature case keeping the Fermi functions instead of 
the step functions. The derivatives with respect to tem- 
perature are calculated from a Sommerfeld expansion. 
An external magnetic field, on the other hand, lifts the 
the degeneracy of the f-level according to 36 
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e f — ► ef + ah . (49) expressed in terms of three spin-dependent integrals 

I oi (uj;a) which which closely parallel their counterparts 
The NCA contribution T,q(u>) — ^ In [e/ + ah — u>] in the absence of an external magnetic field Iq^uj). Keep- 
is directly read off. The Coulomb contribution is now ing only the charge fluctuation contribution yields 



Y,f ] (u) = ~^ UA2 ( -/ oi( w ; cr ) ln I 6 / ~<rh-u>\ + J 2(w; + M^; <r)) 



'01c 



/02c(w;cr) 



f dxf 

o Jo 



(-w+a;+2/)(-w+e/+ ah + x) 



In 



-w + e/— ah + x + y 
>+ef - ah + x + y+ 1 



-f02c(w; Cr) = Io2c(w; 0")/(w; -O") 



dan 



/o {-uJ+ef +ah +x) 

I 03c (uj;a) = / dx dy dz- 

Jo Jo Jo \~ u 



ln 



-uj + 1 + x 

— w + a; 

21n 



x+y+z+l 
x+y+z 



+x+y)(-uj+ef+ah+x)(-u)+ef-ah+x+y+z) 



(50) 
(51) 

(52) 
(53) 
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FIG. 6. Reduced thermodynamical coefficients from Eq. 48 
as functions of A for tf = —2/3 (diamonds) and e/ = —0.7 
(circles). (a):7^, (b):*^, (c):^^, (d):^ 1 ) 

At low temperatures, we find a finite temperature- 
independent Pauli-like spin susceptibility and a lin- 



ear specific heat indicating a nonmagnetic Fermi liquid 
ground state. The results are displayed in Figure 6(a) 
and (b). The coefficient 7/ of the linear specific heat is 
reduced by the conduction electron interactions reflecting 
the enhancement of the Kondo temperature. The scaling 
7/ ~ 1/Tff can be seen from from Figure 6(c) 38 . Simi- 
larly, the magnetic susceptibility is reduced by the con- 
duction electron interactions. The reduction determined 
here is comparable to the value obtained by Hofstetter 
et al. . Its actual values, however, exhibit deviations 
from universal scaling with the inverse Kondo tempera- 
ture 1 /Tr reflecting the importance of quasiparticle in- 
teractions. This is to be expected from the explicit ex- 
pression for the spin susceptibility calculated to leading 
order in the conduction electron interaction 



Xs{U)lxf ] * 1 - 



d_ 

duo 



log 



OI1 2 





- <) 1 








+ 


( 


1 ^0 

. du J 





9*1% 
dh 2 
d 2 Z ( ° 
L dh 2 



ST, 



K 



,h=0 



(54) 



„(°) 



,h = 0y 



The interaction correction to the spin susceptibility 
consists of two terms where the first one is proportional 
to the change in the Kondo temperature. As its coeffi- 
cient varies proportional to Tr-(O) -2 we expect this con- 
tribution to dominate in the close to integer valence limit 
where Tk(0) becomes small. In this limit we should re- 
cover the typical Kondo scaling. The second term which 
results from the variation with h of the interaction cor- 
rections tends to further reduce the spin susceptibility. 
Within the spirit of Landau's Fermi liquid theory it gives 
rise to a positive Landau parameter F^. 

The deviation from simple scaling is also reflected 
in the Sommcrfcld- Wilson ratio which is reduced by 
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the conduction electron Coulomb interaction (see Figure 
6(d)). Since 



(55) 



this quantity allows us to estimate the Landau parameter 



susceptibility Xs is finite at T = 0. Even in the lowest or- 
der in the inverse degeneracy, the low-temperature prop- 
erties are not universal in the sense that their variation 
with U cannot be accounted for by properly adjusting 
the Kondo temperature. 



VI. DISCUSSIONS 




0.25 



FIG. 7. Reduced thermodinamical coefficients from Eq. 48 
as functions of A for tf = —2/3 (diamonds) and ey = —0.7 
(circles), (a,):nf\ (b):^ 

The deviations from the universal scaling with Tk are 
most strongly evident in the f- valence, nj, and in the 
f-charge susceptibility, Xc displayed in Figure 7. The f- 
valence is slightly decreased by the conduction electron 
Coulomb repulsion as can be seen from the coefficient 
n-p displayed in Figure 7, (a). This behavior reflects two 
competing effects, i. e., (a) the increase in characteris- 
tic energy as demonstrated in Figure 5 which is partially 
compensated by (b) the enhancement of the effective hy- 
bridization. 

The f-charge susceptibility Xc is affected more dramati- 
cally as suggested by the coefficient Xc in Figure 7, (b). 
In the Kondo regime, it is enhanced by the Coulomb 
interaction of the conduction electrons. This enhance- 
ment, however, decreases as we approach the mixed- 
valent regime. In this parameter regime, however, the 
adopted approximation (charge fluctuation contribution 
to selfenergy only) ceases to be valid. 

To summarize, the low-temperature properties of a 
dilute magnetic alloy are significantly modified by the 
Coulomb repulsion of the conduction electrons. The lo- 
cal Fermi liquid properties are preserved in the sense that 
the magnetic contribution to the specific heat varies lin- 
early with temperature, Cj ~ T, and that the magnetic 



We formulated here the NCA equations for an Ander- 
son impurity coupled to interacting conduction electrons. 
Due to the CEI an energy-dependent effective hybridiza- 
tion vertex appears in the usual NCA equations. The 
case of weak CEI was investigated in detail here includ- 
ing calculations of measurable thermodynamical proper- 
ties of a dilute magnetic impurities system. The results 
for a rigid, i. e., energy-independent vertex can serve as 
a guide-line for the case of a magnetic impurity embed- 
ded in a Fermi liquid where the central quantity is the 
quasiparticle t-matrix 39 . 

For weak electron-electron interaction the increase of 
Tk may be understood as resulting from the reduced 
probability of finding doubly occupied and empty lattice 
sites in the correlated conduction electron systems. The 
increased number of uncompensated conduction electron 
spins finally leads to the enhancement of the effective 
hybridization coupling. The analogy between the Kondo 
spin model and the Anderson impurity model in its local 
moment regime is not complete in the case of correlated 
conduction electrons ( see also 16 ). In the former case 
Tk will increase monotonously with U because of the 
enhancement of the exchange interaction. In the latter 
case the process is two-staged, it involves the formation 
of a local moment and its interaction with the conduc- 
tion electrons. As we show here CEI influence thermo- 
dynamical quantities in a non-trivial fashion: the specific 
heat coefficient "/(U) scales with Tk(U) like in the non- 
interacting case while the magnetic susceptibility \ a {U) 
does not. That may be understood because the latter 
quantity depends not only on the low energy excitations 
spectrum as the former but also on the matrix element 
which is influenced by CEI. 

There is some controversy about the scaling with 
Tk(U) of Xs(U). In a somewhat artificial model of 15 
the scaling of Xs{U) with Tk(U) is preserved for not too 
large values of U. In the paper 13 this question is not dis- 
cussed explicitly but it may be judged from the relevant 
plots in 13 that their results do not show the usual scaling 
behaviour of Xs(U) but rather resemble our results. 

The present results are based on the separation of en- 
ergy scales. In so far as the vertex correction does not 
produce a new low frequency scale or does not domi- 
natly contribute to the empry state ( empty-f-state ) 
self-energy we anticipate no qualitative changes when us- 
ing more sophisticated approximations for the vertices 
r^g? (1,2; 3, 4), (Figure 1 c) appropriate for the strong 
correlation regime. We may expect that for sufficiently 
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large U the virtual transition from the f-state to the con- 
duction state will cost too much energy inhibiting the 
Tk -increase and leading eventually to the change in the 
trend 12 ' 13 . For a quantitative treatment of this problem 
there are two visible approaches. One is to introduce 
summation of the infinite subseries of the ring, the lad- 
der and the particle-particle type starting with the 2nd 
order self-energies 31 . The other way is to approximate 
the vertex correction by response functions ( dynamic 
susceptibilities ). The other visible development of the 
theory of the Kondo effect for electronic correlated sys- 
tem in general and high-T c cuprates in particular is the 
use of the antifcrromagnon dynamic susceptibility 40 or 
other phenomenological models which develop a short- 
range order with the virtual breaking of singlet-triplet 
degeneracy in the conduction-electron system ( see also 
Ref. 12 ). In addition, the influence of conduction elec- 
tron interactions on the spectral properties of magnetic 
impurities and their dependence upon the doping are also 
interesting topics for future investigations. 

In conclusion, the NCA theory of an Anderson im- 
purity embedded in a metall with correlated conduction 
electrons is developed and general NCA equations for 
the interacting conduction electrons are obtained. It is 
shown that due to the renormalisation of the hybridiza- 
tion interaction the characteristic energy is increased by 
the weak interactions. The influence of weak conduc- 
tion electron interations on thermodynamic properties of 
magnetic impurities is discussed. 
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APPENDIX A: FOURTH ORDER TERM 



We estimate the variation with the orbital degeneracy 
N of the vertex-corrected boson selfenergy (Figure 1(c) 
assuming a simple band structure for the metallic host. 
The orbitally degenerate conduction bands are modelled 
by tight-binding s-bands the dispersion being determined 
by hoping between nearest neighbor sites on a simple cu- 
bic lattice. The magnetic impurity sitting at the origin 
is surrounded by six nearest neighbors. The wave vec- 
tor dependence of the coupling between the conduction 
states and the strongly correlated local orbital (j, m) is 
given by the canonical structure constant Soo-,jm(k) 41 ' 42 
of the Atomic Sphere Approximation (ASA) to the Lin- 
ear Muffin Tin Orbital (LMTO) method. Here j and m 
refer to the azimuthal and magnetic quantum number of 
the impurity orbital under consideration. 



Let us neglect spin-orbit interaction for a first qual- 
itative discussion. The hybridization matrix element 
V ma (k) does not depend upon a. It is given by 




e+i 



(Al) 



where the argument Rj of the spherical harmonic Y £m de- 
notes the unit vector pointing from the impurity to the 
nearest neighbor sites Rj. The overall length scale S is 
usually chosen as the average atomic radius while Vb(e^) 
is an energy-dependent real prefactor. 

Starting from this form of the hybridization we shall 
first derive an expression for the second order term which 
is subsequently compared to the fourth order average. 

The sum 



(A2) 



is easily evaluated using the addition theorem for spher- 
ical harmonics 



(Yim(R 3 )) Y £m (R f ) 



21+1 

47T 



^(hi-Ri) (A3) 



where N = 21 + 1 is the degeneracy of the impurity level. 
For a simple cubic lattice with 

R J-Rj' = S R j R j ,- 5 R j ,-R j , ( A4 ) 

we obtain for i = 3 

£ V^{k)V ma {k')V {e % )V {e p )N x 

m 

^2 (e-ik-Rj e ik'-Ri _ e -ik-Rj e -ik'-Rj^ ( A5 ) 

As expected, the second order term entering the non- 
crossing diagram 

(J2\v mrT (k)\\ = \V (e)\ 2 N(Y^ M . (l-cos2A-4-)) e 



\V (e)\ 2 N2(J2 M . sin 2 k-Rj) e 
\V (e)\ 2 2zN(sm 2 k x ) e 



(A6) 



is proportional to the degeneracy N — 21 + 1 and the 
number of nearest neighbbors with z = 6 for s.c.L. The 
average over the constant energy surface which has to 
be numerically evaluated is of order unity and varies 
smoothly with the energy e. To summarize: The non- 
crossing diagrams involve the combination 



NA(e) = NN{e) \V (e)\ 2 2z(sin 2 k x ) e 



(A7) 
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The contributions from the Coulomb-corrected vertex, 
on the other hand, require the fourth-order term from 
Eq. (11): 

m fel,fe 2 ,fe3,fc4 

E ^k 2 - ^)<% 4 - E 4 ) X 

m,m' m' 

<5*(fci+fc 2 -fc 3 -fc4) (A8) 

where S*(ki + k 2 — k 3 — k 4 ) is the Laue function, Eq. 
(12). If we were to neglect momentum conservation this 
expression would vanish identically. 

Inserting the explicit expressions for the hybridization 
matrix elements Eq. (Al) reduces to a sum of local con- 
tributions 

]T V#> {Ex ,E 2 ,E 3 ,E 4 )=J2Y. ^ E ^ x 

m R n m > m ' 

I m >a 2 {Rn, E 2 )I* na;s (Rn,E 3 )I* ll(J4 (R n , E 4 ) (A9) 

where the averages 

I ma {Rn,E) = N(E)(V ma (k)e^ R -) E (A10) 

are given by 

I m ,(R n ,E) = N(E)V (E) H) V^Y^Y 3 * m (4-) 

( e ik-Rj e ik-ii n j E ( AU ) 

In the next step, we sum over the magnetic quantum 
numbers m and m! 



^ C, (-Rn, E^I^jRn, E 3 ) 

n=-£ 

= V {E 1 )V {E 3 ) (e rk ^ + ^)) El (e-^i^' + ^)) E3 



Rj , R' 



x 4tt ^ y; m (4) r 3 * m fa) 

m=-l 

= JVV (£i)Vo(^) E (e^-^ + ^)) Sl 



(A12) 



The local term Rj vanishes identically due to the sym- 
metry of the constant energy surface. For finite R n , the 
averages over the constant energy surfaces decay with 
increasing R n due to the oscillatory behavior of the in- 
tegrand. The leading contribution to the lattice sum in- 
volves the nearest neighbor term, e. g. R n = (1,0,0). 



In the subsequent summation over the nearest neigh- 
bors Rj, the nonvanishing contributions are given by 
Rj = (-1,0,0) and Rj = (1,0,0) yielding 

(1 - cos 2k x ) El (1 - cos 2k x ) E2 = 4(sin 2 k x ) El (sin 2 k x ) E2 

(A13) 

The leading contribution to the fourth-order term 
hence factorizes according to 



1-J2V^(E 1 ,E 2 ,E 3 ,E 4 ) = 

m 

4 

zN Y[ N(E r )V (E r )2(sin 2 k x ) Er 



(A14) 



r=l 



The average 2 (sin 2 k x ) E is of order unity as mentioned 
above. If we neglect the dependence of the density of 
states and Vo on the energy E we find for the averaged 
value of VW(E U E 2 ,E 3 ,E 4 ) 

,«> = ^(^) 2 (A 15 ) 

In conclusion, we showed that in the leading approx- 
imation of the hybridisation matrix element expansion, 

Eq. (Al), the ratio j^WJJJJJTp — ' lmQm °^ *he or d er 01 
unity. 
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